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Abstract 

In this article we present a systematic derivation of the Maxwell-Bloch equations describing amplification and laser action 
in a ring cavity. We derive the Maxwell-Bloch equations for a two-level medium and discuss their applicability to standard 
three- and four-level systems. After discusing amplification, we consider lasing and pay special attention to the obtention of 
the laser equations in the uniform field approximation. Finally, the connection of the laser equations with the Lorenz model is 
considered. 



I. INTRODUCTION 

Laser theory is a major branch of quantum optics 
and there are many textbooks devoted to that mat- 
ter or thatpay a special attention to it (see, e.g., 

HEHSIalEHIiHOiJIIIl E3) In s P ite ° f this 

we do think that there is room for new didactic presen- 
tations of the basic semiclassical laser theory equations 
as some aspects are not properly covered in the standard 
didactic material or are scattered in specialized sources. 
The most clear example concerns the uniform field limit 
approximation |13| . which is usually assumed ab initio 
without discussion, and when discussed, as e.g. in 0, 
it is done in a way that admits relevant simplifications. 
In fact, this important approximation has found a cor- 
rect form only recently [l4[ . Other important aspect that 
is usually missed in textbooks is the applicability of the 
standard two-level approximation to the more realistic 
three- and four-level schemes. Certainly this matter is 
discussed in some detail in ^1 but we find it important 
to insist on this as it is usually missed and may lead to 
some misconceptions, as we discuss below. 

There are many good general textbooks on the funda- 
mentals of lasers, e.g. 0, [HIS El > an d we refer the reader 
to any of them for getting an overview on the general 
characteristics of the different laser types. Here it will 
suffice to say a few words on the structure of the laser. A 
typical laser consists of three basic elements: An optical 
cavity, an amplifying medium, and a pumping mecha- 
nism, see Fig. 1. The optical cavity (also named res- 
onator or oscillator) consists of two or more mirrors that 
force light to propagate in a closed circuit, further im- 
posing it a certain modal structure. There are two basic 
types of optical cavities, namely ring and linear, that dif- 
fer in the boundary condition that the cavity mirrors im- 
pose to the intracavity field. In ring resonators the field 
inside the cavity can be described as a traveling wave . 
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FIG. 1: Scheme of a typical laser in a Fabry-Perot cavity. 



Contrarily, in linear (also named Fabry-Perot-type) res- 
onators, the field is better described as a standing wave, 
which requires a more complicated mathematical descrip- 
tion than the traveling wave case. 

The amplifying medium can be solid, liquid, gas, or 
plasma. Nevertheless, most cases are well described by 
considering that the amplifying medium consists of a 
number of atoms, ions or molecules of which a number of 
states (energy levels), with suitable relaxation rates and 
dipolar momenta, are involved in the interaction with 
the electromagnetic field. It is customary to adopt the 
so-called two-level approximation, i.e., to assume that 
only two energy levels of the amplifying medium are rel- 
evant for the interaction. Actually a minimum of three 



1 The traveling wave will propagate with a given sense of rota- 
tion, say, clockwise. However, one could in principle expect a 



second field propagating counter— clockwise, i.e., one could ex- 
pect bidirectional emission in a ring laser. Nevertheless, this is 
usually (although not always) avoided by using some intracavity 
elements such as, e.g., Faraday isolators. In any case, we shall 
not deal here with bidirectional emission. 
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or four levels are necessary in order to obtain population 
inversion, and we discuss below how the two-level theory 
applies to these more complicated level schemes. 

Then there is the pumping mechanism. This is highly 
specific for each laser type but it has always the same 
purpose: Creating enough population inversion for laser 
action. When modeling radiation-matter interaction in- 
side the laser cavity one can usually forget the specifics 
of the pumping mechanism (whether it is an electric cur- 
rent or a broadband optical discharge or whatever) and 
describe it through a suitable pumping parameter. In 
this point, the consideration of two-, three- or four-level 
atomic schemes turns out to be important, as it is here 
where the pumping mechanism affects the mathematical 
description as we show below. 

Laser physics studies all of these aspects of lasers but 
here we shall not deal but with the mathematical descrip- 
tion of the interaction between light and matter inside the 
laser cavity. In this article we shall provide a systematic 
derivation of the semiclassical laser equations for an im- 
portant and simple case: The homogeneously broadened 
ring laser, which plays the role of a paradigm in laser 
physics. We shall not consider the important issues of 
inhomogeneous broadening or linear optical resonators, 
because we want to maintain the derivation as simple as 
possible (but not more!). 

As stated, we shall use semiclassical theory, i.e., we 
shall consider a classical electromagnetic field in interac- 
tion with a quantized medium. The quantization of the 
medium is necessary in order to correctly describe ab- 
sorption and amplification as the classical theory (that 
models matter as a collection of forced and damped har- 
monic oscillators) cannot be used for that. With respect 
to the quantization of the field, it is not necessary if one 
(i) is not interested in the field fluctuations, and (ii) ac- 
cepts a heuristic description of relaxation phenomena (in 
particular of spontaneous emission). In any case, the 
quantum theory of the laser requires the use of com- 
plicated mathematical techniques and falls outside the 
scope of our interests here. 

After this introduction the rest of the article is orga- 
nized as follows: In Section II we derive the field equa- 
tion; in Section III we derive the matter equations for 
two-, three-, and four-level atoms or molecules; and in 
Section IV we connect these with the field equation and 
write down the Maxwell-Bloch equations. Then Sections 
V and VI are devoted to the analysis of amplification and 
lasing, respectively. In Section VII we present a clear 
derivation of the uniform field equations, and in Section 
VIII we present the "Lorenz" from of the laser equations. 
Finally, in Section IX we present our conclusions. 

II. THE FIELD EQUATION 

Maxwell's equations for a nonmagnetic material with- 
out free charges yield the wave equation 

V 2 E-c" 2 <9 2 E + V(V-E) =Mo<9 t 2 P. (1) 



Along this article we shall assume that the electric field 
E is a plane wave propagating along the z axis, and write 
it in the form 

E (r, t) = \e£ (z, t) e <(*»-<"*) + c .c, (2) 

where e is the unit polarization vector (fixed polarization 
is assumed), and 

k = lu/c. (3) 

We note that ui is an arbitrary reference (carrier) fre- 
quency. For instance, if light is perfectly monochromatic 
of frequency w we can still choose w ^ w as we allow 
the complex amplitude £ (z, t) to be time and space de- 
pendent. The situation is even clearer when dealing with 
light whose spectrum has some finite width: In this case 
even the concept of "light frequency" is ill-defined, and 
clearly oj can be chosen arbitrarily. 

Given the form for the electric field, by consistency 
with the wave equation, the polarization P must read 

P (r, t) = \eV (z, t) e l ( fcz -^) + c .c. (4) 

Now one must substitute these expressions into the 
wave equation and perform the Slowly Varying Envelope 
Approximation (SVEA) that consists in assuming that 

d 2 U < Lud t U < w 2 C7, (5) 
d 2 U < kd z U < k 2 U, (6) 

for U = £ or V . The physical meaning of this impor- 
tant approximation is clear: One considers that temporal 
(spatial) variations of the amplitudes U contain temporal 
(spatial) frequencies that are much smaller than the car- 
rier frequency (wavenumber) . In other words: The ampli- 
tudes U are assumed to vary on time (space) scales much 
slower (longer) than the optical frequency (wavelength). 
Obviously this approximation excludes the (limit) case 
of ultrashort pulses containing only a few cycles of the 
field, but it is overall very accurate in general, even for 
short pulses as soon as a sufficient number of cycles enter 
within the pulse width. 

After performing the SVEA and multiplying the re- 
sulting equation by c 2 /2iuj, one readily obtains 

(d t + cd z )£ = i^V, (7) 
2e 

which is the field equation of interest. Let us remark 
that the SVEA is a fundamental approximation in laser 
theory as it allows to transform the original wave equa- 
tion, which is a second-order partial differential equation 
(PDE), into a first-order PDE. 

Now we need to calculate the source term V and we 
do this in the next section. 



III. THE MATTER EQUATIONS: OPTICAL 
BLOCH EQUATIONS 

The wave equation (7J relates the slowly varying elec- 
tric field amplitude £ with its source, the slowly varying 
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FIG. 2: Schematic of the two-level atom energy levels includ- 
ing relaxation rates (see text) 



polarization amplitude V . We discuss in this section how 
this last quantity is determined. First we introduce the 
two-level atom model and derive the evolution equation 
for its density matrix, the so-called optical Bloch equa- 
tions. Next, the density matrix is shown to yield the 
information necessary for computing V what allows to 
write a closed set of equations describing the coupled evo- 
lution of field and matter, the Maxwell-Bloch equations. 
Then we consider the case of three- and four-level atoms, 
which is a more realistic approximation to actual lasers. 
After deriving their corresponding Bloch equations, we 
discuss the conditions under which the two-level model 
can be applied to three- and four-level atoms. In partic- 
ular this is a necessary step for correctly understand the 
meaning of the pump parameter. 



A. The two— level atom model 

1. Hamiltonian 

The Hamiltonian of the system consists of two pieces: 
One describing the atom or molecule in the absence of 
electromagnetic interaction, and the other describing the 
action of the electromagnetic field on this atom, i.e. 



JT (r,t) = flat + flint M)- 



(8) 



The material medium is assumed to be a system of iden- 
tical two-level atoms or molecules, i.e., it is assumed that 
the material medium is homogeneously broadened. We 
denote by |1) and |2) the lower and higher energy levels, 
respectively, and by w%i the transition frequency of one 
of these atoms. This means that the atomic Hamiltonian 
flat verifies 



flat |2) = +ifiW2i |2) , 

fl at |l) = -\fvjj 21 \l) , 



(9) 
(10) 



where we have chosen the arbitrary (and unimportant) 
energy origin in such a way that it lies halfway between 



both states energies (see, Fig. 2). The matrix represen- 
tation for this Hamiltonian thus reads 



flat. = 











(11) 



where the level ordering has been chosen to be {|2) , |1)}. 

The interaction Hamiltonian H mt is taken in the elec- 
tric dipole approximation. Roughly speaking, this ap- 
proximation is valid when the light wavelength is much 
longer than the typical dimensions of the electronic cloud, 
which is on the order of 1 A. Thus the approximation is 
justified in the infrared and visible parts of the spectrum 
and even in the ultraviolet. This interaction hamiltonian 
reads 



flint (M) = -A-E(r,t), 



(12) 



where r denotes the position of the atom (which is not 
quantized in the theory) and the operator p, — — ef at , 
being — e the electron charge and f at the vector position 
operator of the electron relative to the point-like nucleus. 
p, acts on the atomic variables whereas in this semiclas- 
sical formalism the field E is a c-number. In the chosen 
basis ordering the matrix form for this hamiltonian reads 



flnt M) - 



"A i 22 
Ml2 



EM) 
EM) 



-M21 

"Mil 



E(r,t) 
EM) 



(13) 



where the matrix elements 



(m| fi \n) 



d 3 r at lj)* m (f at ) fatV'n (?at) , 



(14) 

and ip n (f a t) is the wavefunction (in position representa- 
tion) of the atomic state |n). (Note that jtx m „ = M^ m .) 
We now recall the parity property of atomic eigenstates: 
All atomic eigenstates have well defined parity (even or 
odd) due to the central character of the atomic poten- 
tial. This means that = /x 22 — and then, in order 
to have interaction, we must consider states |1) and |2) 
with opposite parity (this is the basic selection rule of 
atomic transitions in the electric dipole approximation). 
Hence the interaction hamiltonian (|13|) becomes 



flnt M) 



V(r,t) 
V*(r,t) 



where we have introduced the notation 
V(r,t) = -/i 21 -EM) 



(15) 



(16) 



Taking into account the form of the electric field, Eq. J2J, 
V (r, t) becomes 

V (r, t) = ~ha (z, t) e «(**-<"*) - hf3 (z, t) e -*(**-<"*) , (17) 

where we have defined 



a (z, t) 



M21 



2h 



-E(z,t), 0(z,t) 



M21 



2h 



-£*(z,t). 
(18) 
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We note that 2a is usually referred to as the (complex) 
Rabi frequency of the light field. 

Finally the total hamiltonian H, Eq. JHJl, for a two- 
level atom located at position r interacting with a light 
field reads 



H(r,t) 



V* (r,t) 



V(r,t) 



(19) 



2. The density matrix. Evolution 

The hamiltonian H can serve us to write down the 
Schrodinger equation for the atomic wavefunction. In- 
stead, we use here the density matrix formalism as it is 
the most appropriate in order to incorporate damping 
and pumping terms into the equations of motion, some- 
thing we shall do in the next subsection. In the chosen 
basis ordering, the density matrix p representing a two- 
level atom located at r takes the form 



P22 (r, t) 
P12 (r,t) 



P21 (r,t) 
Pn (r,t) 



(20) 



The meaning of the matrix elements is as follows: p m m 
denotes the probability (0 < p mm < 1) that the atom 
occupies state |m), and p mn (— p* im ) is the coherence 
between the two atomic states, which is related with 
the polarization induced in the the atom by the light 
field; see below. The evolution of p is governed by the 
Schrodinger-von Neumann equation 



ihd t p= [H,p] 



(21) 



Upon substituting Eqs. lj2"U|) and into Eq. J3TJ one 
obtains a set of equations which is simplified by defining 
the new variables 



012 = cr 21 = pi 2 e 



i(kz-ujt) 



(22) 



This is motivated by the functional dependence of the 
nondiagonal elements p±2 and p 2 \ on space and time un- 
der free evolution (V — 0). (We note that the above 
transformation is equivalent to working in the so-called 
interaction picture of quantum mechanics.) The explicit 
space-time dependence added in Eq. 12 I'll makes that 
the new quantities cr^ are slowly varying as will become 
evident later. In terms of these reduced density matrix 
elements, and making use of Eq. 1|17|) . the Schrodinger- 
von Neumann equation (|21J) becomes 



dtP22 = iao\2 + i(3<Ti2e 



-2i(kz—ut) 



+ C.C., 



(23) 



dtp n = ~iaa 12 - i(3a 12 e^ kz ~^ + c.c, (24) 
d t a 12 = -iSa 12 + i (p 22 - Pll ) [a* + /3* e ^(^-*)] , 

(25) 

where we have introduced the mistuning, or detuning, 
parameter 



Note that dt (p 22 + Pn) — 0, what implies probability 
conservation. 

We now make a most important and widely used 
approximation in quantum optics, namely the Rotat- 
ing Wave Approximation (RWA). An inspection of Eqs. 
(|23|I - H25|) shows that, in the absence of interaction (£ = 
0, i.e. a — (3 — in the new notation), p 22 and p\\ are 
constant and er 12 = cr^i oscillate at the low (non optical) 
frequency 6. This means that the time scales of the free 
system are large as compared with the optical periods. 
Now, if the interaction is turned on we see in Eqs. I|23|) - 
(12511 that slowly varying terms (those proportional to a 
or a*) appear, as well as high frequency terms oscillat- 
ing as exp [±2i (kz — tot)] (the terms proportional to (3 or 
f3*). Clearly the atom cannot respond to the latter and 
one can discard them. This is the RWA, which can be 
easily demonstrated by using perturbation theory. 

After performing the RWA, Eqs. (|25|> - l|2"5Tl become 

dtP22 = i {acr 12 - a*a 21 ) , (27) 
QtPw = -i (ao-12 - a*a 2 i) , (28) 
d t o-\ 2 = -iSau + ia* (p 22 - pn) , (29) 

which is the standard form of the optical Bloch equations 
for a single atom. 



3. The population matrix 

We are dealing with a situation in which there is not 
a single atom or molecule interacting with the light field 
but a very large number of them, and then some ensem- 
ble averaging must be performed. The ensemble aver- 
aged density matrix is called the population matrix 
although the name density matrix is more frequently 
used obscuring the differences between the two opera- 
tors. Here we are not going to introduce the population 
matrix rigorously and we refer the interested reader to 
[l| or |8j for further details. 

The population matrix of an ensemble of molecules is 
defined as 



p{z,t)=N- l Y, a Pa{z,t). 



(30) 



Here 2 p is the population matrix, p a is the density matrix 
for an atom labeled by a, and a runs along all molecules 
that, at time t, are within z and z + dz . J\f is the number 
of such molecules, which is assumed to be independent 
of z and t. The equation of evolution of the population 
matrix has two contributions: One of them is formally 
like the Schrodinger-von Neumann equation governing 



8 = U> — U>21- 



(26) 



2 Here we are considering a plane wave laser beam propagating 
along z; hence atoms are grouped according to that coordi- 
nate. In the general three— dimensional case, a population matrix 
p (r, t) must be defined at every differential volume, analogously 

to onj. 
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the evolution of the density matrix of a single atom, and 
the other one describes incoherent processes (i.e. not 
due to the interaction with the electromagnetic field such 
as pumping and relaxation phenomena due to collisions 
between atoms or to spontaneous emission) 0] 

d t pi j = (ih)- 1 [H,p] ij +(tp) (31) 

j = 1,2). In Eq. 1|31|) the term Tp is the one de- 
scribing incoherent processes and T is the Liouville (su- 
per)operator. 

Consider the situation depicted in Fig. 2. It corre- 
sponds to the following matrix elements for the operator 

(ty 

{^P) 22 = -72P22 + 712/011 + A2, 

(fp) 21 -(rp-)>-7,P2i, (32) 

(f p) ^ = -71P11 + 721P22 + Ai , 

where 

72 = 7^ + 72!, (33) 

7i=7f 4 + 7i2. (34) 

In the above expressions, 7y describes the relaxation rate 
from level \i) to level \j) (that is, the pass of population 
from level \i) to level \j) due to collisions), and jf xt the 
relaxation rate from level \i) to some other external level 
(see Fig. 2). The term is the pumping rate of level \i), 
i.e., it describes the increase of population of level |z) due 
to pumping processes. Notice that it is not specified from 
where this population is coming as only the dynamics of 
the two lasing level populations is being described. We 
shall come back to this important point in the following 
subsection. 

The value of the different decay constants appearing 
in Tp depend strongly on the particular substance and 
operating conditions. In any case it is always verified 
that 

7-L > ^ (72 +7i) . (35) 

what reflects the fact that the coherence p^ is affected 
not only by the relaxation mechanisms affecting the pop- 
ulations, but also by some specific collisions, known as 
dephasing collisions, which do not affect the populations. 

With the above form for the Liouvillian, the population 
matrix equations of evolution read 

dtP22 = -72P22 + 712P11 + A 2 + i (aau - a*a 2 i) , (36) 
dtPn = -71P11 + 721P22 + Ai - i {a<7 12 - a*a 2 i) , (37) 
<9 t cri2 = - (71 + iS) (T12 + ia* (p 22 - pn) , (38) 

where we have removed the overbar in order not to com- 
plicate unnecessarily the notation. Now pa can be un- 
derstood as the fraction of atoms occupying level i.e., 



it is the population of this level. Notice that in Eqs. i|36l 
EEJ> d t (P22 + Pu) 7^ in general, what reflects the fact 
that the system formed by the atomic levels |2) and |1) 
is an open system in which population is gained and lost 
through incoherent processes. 

In the two-level laser model, internal relaxation pro- 
cesses (those governed by 721 and 712) are usually ne- 
glected, and it is further assumed that the two lasing 
levels relax to the external reservoir with the same rate 
7| I = 7f :rt = 7i a:t ■ It i s eas y t° see that in this simplified 
description of relaxation processes, the pumping rates 

P° 

A* = ^, (39) 

7|| 

with p { \ { the population of level |z) in absence of 
fields (a = 0). Moreover, in this particular case 
dt (p22 + P11) = and then a single equation is needed 
for the description of the populations evolution. The 
population difference is then defined 

d = p22~pn, (40) 
and Eqs. (|36I38|) simplify to 

d t d = 7|| (d Q - d) + 2i (aa 12 - a*a 2 \) , (41) 

d t ai2 = - (7_l + iS) <T\2 + ia*d, (42) 

where 

* = A&-Pii, (43) 

is the population difference in the absence of fields, that 
is, the pump parameter. This is the simplest way of 
modeling pumping. Clearly, do > implies an inverted 
medium (with a larger number of excited atoms than of 
atoms in the fundamental state). If pumping is absent 
do = — 1. Notice that do appears as a free parameter, 
that we can take positive or negative, although we have 
not discussed yet how could it be controlled. 

4- Rate Equations 

It is interesting to write down Eqs. I|3t)l38|l when 
7_L 3> 7| 1 , 5 as in this case the adiabatic elimination of the 
atomic polarization is justified (see Appendix 1). This 
adiabatic elimination consists in making dt<7i2 — 0, and 
then Eqs. (|36I38I) reduce to 

dtP22 = A 2 - 72P22 + 712P11 - R (P22 - Pn) , (44) 
9 t pn = A x - 71P11 + 721P22 + R (P22 - Pn) , (45) 

with R = 2|a| 2 /7_L. 

These equations are known as rate equations and are 
widely used in laser physics as in most laser systems the 
condition for adiabatic elimination is met. Let us remark 
that rate equations describe appropriately the interaction 
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FIG. 3: Schematic of a three-level atom. R represents the 
incoherent pumping and the laser field interacts with the 
|2) — > |1) transition. The arrows indicate decay processes. 



of a light field with a two-level system in two limiting 
cases: When the atomic polarization can be adiabatically 
eliminated, as we have discussed, and also when the field 
is broadband (i.e., incoherent) in which case the factor R 
has a different expression to the one we have derived but 
again depends on the square of the field amplitude |15| . 
We shall make use of these equations in the following 
subsection. 



B. Three— level and four— level atom models 

As we already commented in the introduction, actual 
lasers are based on a three-level or four-level scheme 
rather than a two-level one, the extra levels describing 
the reservoirs from which the pump extracts atoms and 
to which damping sends atoms. In fact this extra levels 
are necessary for obtaining population inversion (do > 0) 
which is a necessary condition for amplification and las- 
ing, as we show below. Although these extra levels do 
not participate directly in laser action 3 , the description 
of its indirect participation is essential in order to cor- 
rectly describe pumping and decaying processes. Here 
we shall derive the Bloch equations for three-level and 
four-level atoms interacting with a laser field and an in- 
coherent pump, and connect these equations with the 
two-level laser equations derived in the previous section. 



ions that are the active ions in Ruby and Erbium lasers, 
respectively. On these ions population is excited from 
the lower state |1) to the upper state |3) by the pumping 
mechanism. Then population is transferred from level 
1 3) to the upper lasing level |2) (which is long-lived) by 
relaxation processes, which are extremely fast in these 



We shall model the pumping transition |1) 



1 3) via 



rate equations 4 like Eqs. I|44I45(I . and the interaction of 
the monochromatic field with transition |1) < — > |2) with 
the already derived Bloch equations for a two-level atom. 
As for the relaxation processes, we describe them heuris- 
tically, see Fig. 3. Then we can model these processes 
with the following set of Bloch equations 

9tP33 = ~ (731 + 732) P33 + R [pii ~ P33) , (46) 

dtP22 = -721P22 + 732P33 + i [a<Ti2 - a*a 2 i) , (47) 

dtPll = 721P22 + 731P33 + R (P33 - P11) (48) 

- i (ck7i2 - a*a 2 \) , 

d t ui2 = - (71 + i$) °2i + ia* (P22 - P11) , (49) 

where R is the rate at which ions are pumped by the 
incoherent pump field from level |1) to level |3). Let us 
remark that in writing Eqs. (|46|l and (i), we have 

taken into account all possible transitions due to incoher- 
ent processes with suitable relaxation rates as indicated 
in Fig. 3; and (ii), the incoherent pumping of popu- 
lation from level |1) to level |3) is modelled by the term 
R (pn — P33) appearing in Eqs. (|4*rj|) and (|4*5|l with R pro- 
portional to the pump intensity, i.e., we have described 
the interaction of the pump field with the pumped tran- 
sition by means of rate equations similar to Eqs. i|44|) 
and l|45l) but taking Xi — as all incoherent processes 
have been consistently taken into account. 

Let us further assume that 732 ^ 721, 731, 7_l, R, as it 
occurs in usual three-level lasers. Then we adiabatically 
eliminate the population of level |3). By making dtp33 = 
we get 



P33 



Rpn 

732 



(50) 



This equation shows that 732/233 is a finite quantity; that 
is, P33 is vanishingly small in the limit we are consid- 
ering. Then we can neglect p 33 in Eq. l|4*8|) and put 



1. Bloch equations for three-level atoms 

Consider the three-level atom scheme depicted in Fig. 
3, which can be regarded as an approximate description 
of, e.g., the relevant atomic levels of the Cr 3+ or the Er 3+ 



There is a very important exception. In coherent optically 
pumped lasers the pumping mechanism is a laser field tuned to 
the pumping transition. If the atomic coherences cannot be adia- 
batically eliminated, Raman processes are important and cannot 
be neglected. We shall not consider these lasers here. 



4 For example, in the case of Ruby lasers, pumping comes from 
an incoherent light source, namely a flashlamp. Then the inter- 
action of this incoherent light field with the pumping transition 

1) < ► |3) can be described with the help of rate equations. 

The case of Erbium lasers is different: In this case the pump- 
ing is made with the help of a laser field tuned to the pumping 
transition. In spite of the coherent nature of the pumping field, 
a rate equations description for the pumping transition is also 
well suited in this case, because the adiabatic elimination of the 

atomic coherence of transition |1) < > 3) is fully justified as 

its coherence decay rate is very large as compared with the rest 
of decay rates. Other systems have other pumping mechanisms 
(e.g., the pass of an electrical current through the active medium) 
for which the rate equations description os also appropriate. 
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FIG. 4: Dependence of the effective pump parameter do on 
the normalized actual pump strength R/j for three- (full line) 
and four-level (dashed line) lasers. 7 = 721 for the three-level 
laser and 7 = 720 + 721 for the four-level laser. 



732/°33 = ftpn in Eq. I|47|) . By further noticing that af- 
ter the approximation dt (pn + P22) = 0, we can write 
the simplified model 

dd = ft — 721 - (ft + 721) d + 2i (ao\2 - a*a 2 i) , 

(51) 

dto-\2 = - (71. + id) a 12 + ia*d, (52) 

where d = (P22 — Pn)- These are appropriate Bloch 
equations for most three-level systems. We can now com- 
pare these equations that describe three-level atoms with 
Eqs. (|41(l and (|42J) that describe two-level atoms in a 
simple and usual limit. It is clear that they are isomor- 
phic. Then we can conclude that incoherently pumped 
three-level atoms can be described with the standard 
two-level atom Bloch equations by making the following 
identifications 



7|| — >• ft + 721, 

ft - 721 

do — > • 

ft + 721 



(53) 
(54) 



Notice that (i) the decay rate 71 is pump dependent for 
three-level atoms, and (ii) that the pumping rate do de- 
pends in a nonlinear way with the actual pump parameter 
ft. In Fig. 4 we represent do as a function of the actual 
pumping parameter ft; notice that increasing ft in a fac- 
tor ten, say, does not mean to do so in do. Apart from 
this, we have shown that an incoherently pumped three- 
level medium can be described as a two-level one when 
the adiabatic eliminations we have assumed are justified, 
which is the usual situation. 



2. Bloch equations for four-level atoms 

Consider now the four-level atom scheme shown in Fig. 
5. It can be regarded as an approximate description of, 
e.g., the relevant atomic levels of the Nd 3+ ion that is the 
active ion Nd-YAG or Nd-glass lasers. Assuming, as for 




laser 
field 



FIG. 5: Schematic of a four-level atom. R represents the inco- 
herent pumping and the laser field interacts with the 2) — > |1) 
transition. The arrows indicate decay processes. 



three-level atoms, that the pumping field acting on the 
transition |0) — |3) can be described by rate equations, 
we are left with the following optical Bloch equations 



732) P33 + ft(poo 

732 P33 



dtP33 = - (730 + 731 
dtP22 = - (720 + 72l) P22 

+ i (aa 12 — 0*0-21) , 
dtpn = -710P11 + 721P22 + 731P33 

+ i (aa 12 — a*02i) , 
dtPoo = 710/011 + 720P22 + 730P33 - ft (Poo 
dtO-12 = - (71 + iS) 0-21 + ia* (p 2 2 - pn) ■ 



P33J 



(55) 
(56) 

(57) 



P33) (58) 
(59) 



We can now proceed in a similar way as we did with 
three-level atoms: Let us assume that 732 is much larger 
than any other decay rate and adiabatically eliminate 
P33. Now we get 



P33 



Rp, 



'00 



732 



(60) 



and neglecting the terms 731^33,730^33 and ft/333 we are 
left with 

<9 t p 2 2 = - (720 + 72l)P22 + ftpoO (61) 

+ i (ck7i2 - a*a 2 i) , 

dtpn = -710/011 + 721P22 + i (a<7i2 - a*a 2 i) , (62) 

dtPoo = 710/011 + 720P22 - ftpoo (63) 

d t (Ji2 = - (71 + iS) (721 + ia* (P22 - /Oil) • (64) 

Now we must take into account that the lower lasing level 
|1) usually relaxes very fast towards level |0). This means 
that pn ss and consequently that d = P22 — P11 ~ P22- 
Taking this into account and also that poo + P22 ~ 1 in 
this approximation, we are left with 



dtP22 
9*0-12 



: - (720 - 
i (ao\2 
- - (7-L - 



- 721 + ft) d 

- a*a 2 i) , 

iti) 021 + ia*d 



ft 



(65) 
(66) 
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We see that, after the adiabatic elimination of ^33 and 
Pll, the four-level Bloch equations are isomorphic to Eqs. 
(|41|l and (|42|l that describe two-level atoms. Then we can 
apply the two-level description to a four-level atom by 
making the following identifications 



7| 
d 



720 



- 721 
R 



R, 



720 + 721 + R 



(67) 
(68) 



Again, as it was the case for three-level lasers, the pop- 
ulation decay rate 711 and the pumping rate d of the 
two-level theory must be reinterpreted when applied to 
four-level lasers. 

Once we have shown that the two-level theory of Eqs. 
(|41J) and Ij42(l can be applied to three- and four-level 
lasers by suitably interpreting the parameters 71 and do, 
in the following we shall always refer to the two-level 
model but the reader must keep in mind that the trans- 
formations we have derived must be taken into account 
when applying this theory to three- and four-level lasers. 



IV. THE MAXWELL-BLOCH EQUATIONS 



We finally come back to the wave equation J7J), multiply 
it by (fi 2 i ■ e) /2fi, and make use of Eqs. (|18f) and l)72|) 
for obtaining the final field equation, which we write to- 
gether with the Bloch equations (|41|) - (|42[) for the sake of 
convenience 

(> ' ] ' (73) 



da 

at az 



dt(?i2 = - (71 + iS) C12 + ia*d, (74) 
d t d = 7|| (d — d)+2i (aa i2 - a*a 2 i) , (75) 

where we have introduced the radiation-matter coupling 
constant 

JVuj \n 21 ■ e| 2 



9 



2e h 



(76) 



Note that Eqs. i|73H75[l . form a closed set of equations 
that completely determines, self-consistently, the interac- 
tion between a light field (of amplitude proportional to a, 
see Eq. (^HJ) an d a collection of two-level atoms. This 
set of equations is known as the Maxwell-Bloch equa- 
tions for a two-level system, which can be applied to 
three- and four-level systems by introducing the param- 
eter changes I|53I54[) and (|67I68|I . respectively. 



Once the field equation JJJ) and the optical Bloch equa- 
tions for matter dynamics Eqs. (|41|l and (|42|l have been 
derived, we only need to connect them in order to obtain 
a closed set of equations for the analysis of amplification 
and laser dynamics. 

Under the action of the light field each atom devel- 
ops an electric dipole. As the number of atoms con- 
tained in a small volume (small as compared with the 
light wavelength) is always large, one can assume that at 
each spatial position r there exists a polarization given 
by the quantum-mechanical expectation value of the elec- 
tric dipole moment operator p,. When using the density 
(population) matrix formalism, this expectation value is 
computed as the trace 



P{r,t)=AfTr(p(r,t)fi) 



(69) 



where M denotes the number of atoms per unit volume. 
Making use of Eq. (|2U|) and of the matrix form for the 
dipole moment operator 



M21 
M12 



one has 



P(r,t)=Af[n 12 p2i(r,t)+c.c), (70) 
which, making use of definitions 122II . reads 



P(r,t) = Af V 12 o- 21 (r,t) 



+ ex. 



(71) 



which, compared with Eq. yields 



V. AMPLIFICATION 

The simplest issue that can be studied within the de- 
veloped formalism is the amplification of a monochro- 
matic light beam after traveling some distance along a 
medium. If we identify lo with the actual light frequency, 
then £(z,t) — £ (z), see Eq. ©, which implies that 
a (z, t) — a (z). On the other hand, after a short tran- 
sient (of the order of the inverse of the decay constants), 
the atomic system will have reached a steady configura- 
tion, which is ensured by the presence of damping. Thus, 
after that transient one can ignore the time derivatives in 
the Maxwell-Bloch equations. Solving for the material 
variables (|41|I - H42|I in steady state, one has 



d - 



ll+S 2 



S 2 +jI + 4 7± |a| /7|| 



021, 



do a 



Ha 



(5 2 +71+47_lH /7|| 



(77) 



(78) 



(79) 



where the subscript "s" refers to the steady state. Sub- 
stituting the result into the field equation l|73|l one has 

da dog 7j_ + 
— = ^ a. 

az c £ 2 + 72 +4 7± |a| /7|| 

This equation governs the spatial variation of the field 
amplitude a along the atomic medium. 



Weak field limit 



Before considering the general solution let us concen- 
(72) trate first on the weak field limit, defined as \a\ 2 <§; 
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7±7||/4- In this case the last term of the denominator 
in Eq. 1)79(1 can be ignored and the solution reads 



a(z) = a (0) exp 



a/2 



l + i- 



where 



2d g 

C7J_ ' 



. (80) 



(81) 



and g given by Eq. l|7njl. Parameter a is responsible for 
the attenuation (when a < 0, i.e., when do < 0) or ampli- 
fication (a > 0, i.e., do > 0) of the light along its propaga- 
tion through the material. In case of attenuation, the in- 
verse a -1 is known as penetration depth. In case of ampli- 
fication a receives the name of small-signal gain per unit 
length. (Note that for 6 = 0, \a (z)\ 2 = \a (0)| 2 exp (az).) 

On the other hand the imaginary exponent corre- 
sponds to a correction to the light wavenumber. In fact, 
noticing that a is proportional to the field amplitude £ 
and recalling Eq. J5J one has that the actual wavenumber 
is 

fc cff = k + 8k, (82) 
and consequently the refractive index n = ck e s/to reads 
ck d g (V7-l) 



w 7±w i + (5/j±) 

d N\n 21 -e\ 2 (5/71 



1 + 



2e fr7± 1 + (<S/V) 2 ' 



(83) 



which has the same qualitative behavior as the classi- 
cal expression obtained from the (harmonic oscillator) 
Lorentz model 



B. Strong field limit 

In the opposite limit, namely |a| 2 TJ-T| | / 4, <5 2 , Eq. 
(f79ll becomes 



da 

dz 



l\\dog ( 
4c 



7i 



(84) 



Multiplying this equation by a* and taking the real part 
of the resulting equation one has 



d\a\ 2 _ J\\d g _ 7||7_i_a 



dz 

whose solution reads 

2 



2c 



|a (*)!=!« (0)1' + 



7||7_La 



(85) 



(86) 



Again, amplification requires a > 0, i.e., do > 0. This 
result means that, for strong fields, there exists satura- 
tion: the amplification (whenever do > 0) persists but it 
is linear in the propagation distance, differently from the 
weak field limit, in which amplification occurs exponen- 
tially, see Eq. ifgUl) . 
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FIG. 6: Field intensity during amplification as a function of 
the normalized distance z/z p d, with z v d = a -1 the penetra- 
tion depth. The dashed lines correspond to the weak and 
strong field approximations, Eqs. 18011 and 186H . respectively. 



C. General solution 

In order to consider the general case it is convenient to 
use a polar decomposition for a as 



a(z) 



a(z 







(87) 



Substituting this expression into Eq. (|79f) and separating 
it into its real and imaginary parts one obtains 



d\a\ j±d g 



dz C J2 +7 2 +47 2 | Q |2y 7|| 

d(j) Sdog 1 

dz c S 2 +72 +47J |a| 2 /7||' 

Eq. can be integrated to yield 



(89) 



\a(z)\ 



(7i+^)lnB^ + 2^[| 



|a(0)| 7|| 



a(z)| 2 -|a(0)| 2 



7±dog 



-z, 



(90) 

which does not allow an explicit expression for \a (z)\. In 
any case, Eq. (fSl^l shows that \a (z)^ 1 d\a (z)\ /dz has 
the same sign as do, so that do > implies amplification. 
In Fig. 6 the solution of Eq. (|9"0|) is represented as a 
function of z together with the weak and strong field 
approximations derived above. 

As for Eq. I|89|l , the phase can be determined by notic- 
ing that 



d<j)/dz 



6 1 



from which 



d I o: I d |cv| /dz 7j_ \a\ ' 
7-L I a (0)| 



(91) 



(92) 



Note that, on resonance (5 = uj — uj^i = 0) there is no 
phase variation along the propagation direction (apart 
from the original phase kz). 
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output 




z = 



which, upon using Eq. I|94|l and recalling that k = uj/c 
(this was our choice in writing Eq. J2J0? reads 



£ (0,*) = Tle tkL *£ (L m , t - At) . 



(96) 



Finally, multiplying this equation by (/i 2 i ' e ) /2fi and 
recalling Eq. I|18fl. one has 

a(0,t) = TZe ikL ^a(L m ,t- At). (97) 

We analyze next the monochromatic lasing solution. 



FIG. 7: Scheme of the ring cavity. The active medium is 
placed in the region < z < L m . The black cavity mir- 
rors are perfectly reflecting whilst the grey mirror has a finite 
reflectivity 1Z. The arrows indicate the propagation of the 
intracavity and output fields. 



VI. LASING 

Differently from the previous analysis, in which we as- 
sumed that a given field (whose frequency and initial am- 
plitude are known data) is injected into the entrance face 
of a material, the light field in a laser is not fixed exter- 
nally but is self-consistently generated by the medium, 
through amplification, and must verify the boundary con- 
ditions imposed by the cavity. As the model we have 
developed considers a traveling wave (moving in one di- 
rection) the following analysis only applies to ring lasers 
in which unidirectional operation can take place (in lin- 
ear, i.e. Fabry-Perot, resonators there are two counter- 
propagating waves that form a standing wave, a more 
complicated case that we shall not treat here). 



A. Boundary condition 

We assume that the medium is of length L m and that 
the cavity has a length L c , see Fig. 7. We take z — 
as the entrance face of the amplifying medium. The 
boundary condition imposed by the resonator reads 



E(0,t) = 1ZE (L m , t — At) . 



(93) 



where 1Z represents the (amplitude) reflectivity of the 
mirrors (0 < 1Z 2 < 1 gives the fraction of light power 
that survives after a complete cavity round trip) and 



At 



(94) 



is the time delay taken by the light to travel from the 
exit face of the medium back to its entrance face after 
being reflected by the cavity mirrors. 

Making use of Eq. J3J, and after little algebra, the 
boundary condition (|93|l reads 

£ (0, t)=K exp [i (kL m + uAt)] E (L m , t- At), (95) 



B. Monochromatic (singlemode) emission 

We note that the frequency u> appearing in the field 
expression @ is by now unknown. Under monochro- 
matic operation the laser light has, by definition, a single 
frequency. If we take lu to be the actual lasing mode fre- 
quency, the field amplitude must be then a constant in 
time i.e., a(z,t) — a{z), as in the previous analysis. 
Thus Eq. becomes 



a(0) = TZe ltc ^a(L m ) . 
Using now the polar decomposition l|87|) one has 



\a(0)\ 2 
0(0) 

being m an integer. 



TZ 2 \a(L m )\ 2 
(j) (L m ) + kL c 



2TO7T, 



(98) 

(99) 
(100) 



Determination of the laser intensity 



Let us first analyze the laser intensity |cv| . (In fact 
the laser intensity is proportional to \£ | but remind that 
£ oc a, and then \£\ cx \a\ .) We note that, as we are 
dealing with a field whose amplitude is time independent, 
the analysis of amplification performed in the previous 
section is directly applicable. Making use of Eq. 199(1 . 
Eq. (|90|l becomes, for z — L m , 



l±dog 



1 
2 

,7-L 



7 | +S 2 ) lnTT 



+ 2^(l-U 2 )\a(L n 
7|| 

which, after trivial manipulation yields 
I hi ft 21 



| a (in 



7||7J 



[r — 1 — A 



(101) 



(102) 



4 i -n 2 

where we made ln^" 2 = |ln^ 2 | (remind that TZ 2 < 
1) and we have defined two important parameters, the 
adimensional pump r and the normalized detuning A 
through 



2d gL m _ 

7± c|ln^ 2 | ~ |m7£ 2 |' 

U>21 



S OJ 

A = — = - 

7-1 7-L 



(103) 
(104) 



10 



0.05 

0.04 

' 0.03 


3- 0.02 

0.01 



0.00 



_ 1 1 1 


1 i 1 i y. 


— 


0,5/ ~ 










— -rT"T 


0.99 

" 



0.0 



0.2 



0.4 0.6 

z/L m 



1.0 



FIG. 8: Intracavity field intensity as a function of distance 
for 7|| = 7_l, A = 0, r = 1.5 and the values of 1Z marked in 
the figure. 



We note that the adimensional parameter r is propor- 
tional to the gain properties of the medium and inversely 
proportional to the damping properties of the system. 
In fact, aL m gives the small-signal single-pass gain along 
the amplifying medium (remind that a, Eq. (|81l) . is the 
small-signal gain per unit length). Thus r acts as an 
effective pumping parameter, as will become clear next. 
Equation (|102f> determines the value of the field inten- 
sity at the exit face of the amplifying medium. Clearly, 
in order to be meaningful, \a (L m )\ 2 > 0, what implies 



r > r n 



(105) 



Thus parameter r must exceed a given threshold (the 
lasing threshold ro) in order that the laser emits light. 
This is why r is called the pump parameter (there is a 
minimum pump required for the system starts lasing). 
What we have obtained is the field intensity at the faces 
of the active medium, Eqs. <|99|l and 1|1U2[1 . But it is also 
interesting to analyze how this intensity varies along the 
active medium. Thus, after using Eqs. (|99|l . I|103|) and 
(|104f) . we write down Eq. (|90|l in the form 



1 in, V ; K 2 \a(L n 



\a(z)\ 2 ~TZ 2 \a(L m )\' 



7H7-1 



(106) 



2. Determination of the laser frequency 

Even if it can seem that we know the lasing intensity 
value, the fact is that we still ignore the value of the 
lasing frequency lu and thus the value of A. This problem 
is solved by considering the phase boundary condition 
QlUUfl . First we recall Eq. I|92|). which we write in the 
form 



0(i m )-0(O) 



r Aln 



I a (in 



-±A|lnft 2 | 



(107) 



where Eq. I|99|l has been used in the last equality. Com- 
parison between Eqs. I|100fl an d l|107l) yields 



±A|lnft 2 | = 2?rm- kL c 



(108) 



We now introduce the wavenumber k c and frequency lo c 
of the cavity longitudinal mode closest to the atomic res- 
onance. As we are dealing with a cavity longitudinal 
mode, it must be verified, by definition, that 



k c = 2i:m c /L c , 



— ck c 



(109) 



being m c an integer. Substituting these quantities into 
Eq. (|108[> one gets 



±A \lnTZ 2 \ = {lu c -lo] 



27m, 



(110) 



where n = m — m c is a new integer. We finally recall Eq. 
(|104|l so that Eq. I|110[) yields the following value for the 
laser frequency 

KW 2 1+7±W C . 7-1 2 7TC 

u> n = h n — (HI) 



K + 7_L 

where we have defined 



\IylK 2 
2L C 



(112) 



which is known as the cavity damping rate for reasons 
that will be analyzed in the next section. We note that 
Eq. (|lllf> indicates that there exists a family of solu- 
tions (labeled by the integer n). As we show next all 
these solutions have, in general, different lasing thresh- 
olds. From Eq. (|lllfl the lasing threshold l|105|) can be 
finally determined as 



with \a (i m )| given by Eq. H102JI . This equation can be 
solved numerically and in Fig. 8 we represent its solutions 
for fixed parameters and several values of the reflectivity 
1Z 2 , showing that as 1Z 2 approaches unity the solution 
becomes progressively uniform. This fact suggests that 
for 7?. 2 — > 1, it must be possible to rewrite the laser equa- 
tions in a simpler way as in this limit the steady state is 
independent of z. We shall come back to this point in 
the next section. But first we shall continue analyzing 
the laser steady state. 



^21 



7J 



(113) 



Now, the difference between the cavity and atomic tran- 
sition frequencies is obviously smaller than the free spec- 
tral range, i.e., \lu c — uj2i\ < 2t:c/L c . This makes that 
r on is minimum for n = and, also, that the frequency 
of the amplified mode, ujq, be given by 



w 



KUJ 2 1 + J±UJ C 
K + 7-1 



(114) 
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which is the pulling formula. The result is that the laser 
frequency is a compromise between the cavity and atomic 
transition frequency. Notice that for a "good cavity", 
k < 71, the laser frequency approaches the cavity fre- 
quency, whilst for a "bad cavity", 7j_ <C K, the laser 
frequency approaches that of the atomic transition. This 
is a quite intuitive result indeed. 



3. The resonant case 

Let us analyze the relevant case uj c = W21, correspond- 
ing to a cavity exactly tuned to the atomic resonance. In 
this case the pump must verify 



r > r n 



2ttc 



(k + 7_l) L c 



(115) 



and the lasing mode with lowest threshold is that with 
n — 0, as discussed above. Hence, at resonance, the basic 
lasing solution has a threshold given by r on — 1, and its 
frequency is lu = uj c = u>2i, see Eq. (|111H for n = 0. 

The amplitude of this lasing solution verifies Eq. I|79|) 
with 5 = lu — LU21 = 0: 



da s 
dz 



dog 



1 



7_lc 1 



fCt s 



(116) 



7X7|| 



We note that we have introduced the subscript " s" to em- 
phasize that this amplitude a corresponds to the steady 
lasing solution. 

Finally, the "intensity" of the laser light at the exit 
of the active medium is given by Eq. I|102[l with A = 
S/y± = 0: 



\a s (.Lm)!" 



7||7i 



In ft 2 



4 l-W 2 



(r-1). 



(117) 



We note that there exists no phase variation of the laser 
complex amplitude a s along the medium, see Eq. I|89|) 
with 8 — 0. We shall make use of these expressions in 
the following section. 



VII. THE LASER EQUATIONS IN THE UNI- 
FORM FIELD LIMIT 



In this section we want to find a simpler model that 
allows us to study laser dynamics and instabilities in an 
easy way. The sought model is known as the Lorenz- 
Haken model and can be rigorously derived from the 
Maxwell-Bloch equations (glJ-JEJ an d <C3 m tne so- 
called uniform field limit, which we consider now. This 
limit assumes that the cavity reflectivity is closest to 
unity (1Z 2 — > 1 in all previous expressions). For the 
sake of simplicity |lfj the derivation will be done in 
the resonant case, where the cavity is tuned in such a 
way that one of its longitudinal modes has a frequency 
lu c that matches exactly the atomic resonance frequency 



u>2i- In this case the analysis done in Sec. (jVI B 3|) sug- 
gests to choose the value of the arbitrary frequency to as 
to = lo c = u>2i. (We remind that we can choose freely 
this value. If this election is "wrong" the laser equations 
will yield an electric field amplitude a which contains a 
phase factor of the form exp (— iAut) that will define the 
actual laser frequency.) 

First we recall the Maxwell-Bloch equations (|73H75|) 
for S = 0: 



da da 
dt dz 



dt 



-7-L°21 



iad. 



^ = 7|| (d - d) + 2i (ao 12 ~ a* 'a 2 i) , 



(118) 
(119) 
(120) 



which are to be supplemented by the boundary condition 
(|97|l with k — k c — uj c /c, see Eq. @, so that k = 
2irm c /L c , see Eq. (|109|) . With these assumptions the 
boundary condition (|97|l becomes 



a(0,t) =TZa(L m ,t- At). 



(121) 



We note that this boundary condition is not isochronous 
(it relates values of the field amplitude at different times) 
and this makes difficult the analysis. We note for later 
use that this boundary condition applies, in particular, to 
the steady lasing solution (independent of time) so that 

a s (0) = Ka s (L m ) . (122) 

These equations form the basis of our study. 



A. A first change of variables 

In order to make the boundary condition isochronous, 
we introduce the following change of variables pj: 



a' (z, t) = a(z,t — t) , 

°2i ( z ,t) = °2i (z,t-r) , 
d' (z, t) = d(z,t — t) , 



(123) 
(124) 
(125) 



with r = zAt/L m and At = (L c - L m ) /c, Eq. (jHU- The 
new variables verify 

da da' , . 

at = -W' (126) 

da da' . At da' 

(127) 



dz 



dz 



L m dt 



and similar expressions for the material variables. Sub- 
stitution of the previous relations into Eqs. (|118fl - H12(J|) 
yields 



L, dc 



da' 



dt 



= ~ ia'd' , 



dd' 

-Q^T = ~7|| ( d> - d o) + 2iaa 12 + c.c. 



(128) 
(129) 
(130) 
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where we used Eq. (|94() . According to Eq. I|121|) the new 

variables verify the following boundary condition 



a' (0,t) = Ka' (L m ,t) , 



(131) 



which is now isochronous. We note that the definition of 
the new variables is mathematically equivalent to " bend" 
the active medium on itself so that its entrance (z = 0) 
and exit faces (z = L m ) coincide. 



B. A second change of variables 

Now we define another set of variables by referring the 
previous ones to their monochromatic lasing values ana- 
lyzed in the previous sections. The steady values of the 
material variables have been calculated in Sec. 10, Eqs. 
()77J) and H78J) . that, particularized to the case S = we 
are considering, read 



d s (z) = d Q 



1 



0-21,s [Z) 



7-L7|| 1 S| 

-idoa s 1_ 

7-1 1 



(132) 
(133) 



7-L7| 



We note that these quantities are z— dependent as a s is, 
Eq. Ullfijl . In particular we define the new variables 
through 



„ / s i a' (z, t) 

F(z,t) = V^T^-Y, 
a s (z) 

P( M ) = Vr=T^iM 

0-21,8 (Z) 



D(z,t) 



d'(z,t) 
d s (z) 



(134) 
(135) 
(136) 



where a s verifies Eq. (|116[) . The equations for F,P,D 
are obtained from Eqs. I|128|I - 1|13U[) . First the equation 
for F is computed. From its definition we have 



dF , 1 da' 

= V r ~ 1 ST' 

at a s at 



dF 



1 da' a' 1 da^ 



dz \ a s dz a s a s dz 



(137) 
(138) 



Making use of these and of Eq. 11161) we build the fol- 
lowing equation for F 



L c dF dF Vr~=T . 

L m dt dz a s 



dog 



1 



7± l + ^-b s | 

7X711 1 1 



;F. (139) 



that by using the definition of P and of Eq. (|133fl trans- 
forms into 



dF dF 



(140) 



where 



C F (z) = 



cL m 
d g L r 



7± L c i 



4 I |2 • 
-TT- \aJ 



(141) 
(142) 



(Note that v < c has the dimensions of a velocity.) The 
equations for the material variables are easier to be ob- 
tained. Making use of the definitions of F, P and D, and 
making use of Eqs. (|129|l and (|130J) we obtain 



dP 

~b~t 

dD 
~dt 
where 



= — 7j_P + CpFD, 



(143) 



- ll] (D-D )-(C D FP* + C* D F*P), (144) 



C P {z) 



021,! 



d R 



C D (z) = 



-2ia s a2 1 . 



(145) 
(146) 
(147) 



d s (r - 1) ' 

After using the steady state equations i|132|) and 11331) 

these expressions can be written as 

Cp = 7i, 



D (z) = l-\ \a s \ , C D (z) = — . 

7-L7|| 7±( r -l) 



(148) 
(149) 



Up to this point, the equations for F, P, and D are 
equivalent to the original Maxwell-Bloch equations, as 
no approximation has been done. 



C. The Uniform Field Limit 

We now study the behavior of Cf, Cd and Do in the 
case when the cavity mirrors have a very good qual- 
ity, i.e., when the reflectivity 1Z is very close to unity. 
In this limit, the boundary condition (|122f) says that 
a s (0) « a s (L m ). On the other hand, the steady state 
equation l|116|l tells us that |a s | is a monotonic grow- 
ing function of z. Under these circumstances, one can 
assume, to a very good approximation, that \a s (z)\ 2 is 
a constant along the amplifying medium. (These facts 
can in fact be seen in Fig. 8.) In this case its value co- 
incides, for instance, with its value at the medium exit 
face, |a s (L m )| 2 , which is given by Eq. (|117fl : 



\a s (z)!' 



7||7i 



In ft 2 



4 1-K 2 



1) 



Vz. 



But, as we are considering the limit TZ — * 1, the quotient 
|ln7£ 2 | / (l — TZ 2 ) also tends to unity as can be checked 
easily, and we finally have 



7||7i. 



(r-1). 



Vz. 



(150) 
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This space uniformity of the laser intensity along the am- 
plifying medium when 72. — >■ 1 is the reason for the name 
" Uniform Field Limit" . (We note that in the literature 
the uniform field limit has been customarily associated 
not only with the high reflectivity condition but also with 
the small gain condition aL nl — * 0. We see here that the 
latter condition is completely superfluous.) Substitution 
of ifTCTjl into Eqs. dUl and lfTl^|) yields: 

d gL m 
C F (z) = — , 

D (z)=r, C D (z) = ^, Vz. (151) 

Finally, making use of definitions (|1U3[) and (|112H . Cp 
simply reads: 

C F (z) = K, Vz. (152) 



we have 

2 

F(z,t) = a'(z.t), 

P(z,t) = 2i^- f^a' 21 (z,t), 
do V 7|| 

D(z,t) = ^-d! (z,t) . 
do 

Thus F has the meaning of a laser field amplitude, P 
has the meaning of material polarization and D has the 
meaning of population difference. 

Equations (|153ll 1155|) allow to study two types of laser 
operation: singlemode and multimode. The method for 
deriving the laser equations in the uniform field limit 
we have followed here was presented in ^4| (see also 
where an application to multimode emission was 
addressed. In the following we shall concentrate on the 
singlemode laser. 



1. The Laser Equations in the Uniform Field Limit 

Substitution of expressions (|148|) , 11511) and (|152[) into 
Eqs. CU, JHI and yields 



dD 
~dt 



dP 

-^=^(FD-P), 
7|| [r- D~\{FP* +F*P)] 



(153) 
(154) 
(155) 



We finally need to consider the boundary condition that 
applies to these equations. By considering the definition 
(Pljl for F 

F(0,*) = ^,F(W) = ^%f 
a s (0) a s (L m ) 

Making the quotient of both quantities we have 

F(0,t) = a s (L m ) q'(0,t) 
F{L m ,t) a s {0) a'{L m ,ty 

which, making use of Eqs. I|122fl and i|131|) yields 



F (0, t) = F (L m , t) . 



(156) 



We thus see that the boundary condition for the field 
amplitude is periodic (we note that this is not due to 
the uniform field limit but to the very definition of F). 
This is of great importance as will allow us, owing to the 
Fourier theorem, to decompose F in terms of periodic 
functions. 

Before studying the obtained Maxwell Bloch equations 
(|153fl 1155fl . let us demonstrate that F, P, and D, are 
equivalent to the original variables, apart from constant 
scale factors. From definitions i|134|) . and using the uni- 
form field limit results developed in the previous section, 



VIII. THE SINGLE-MODE LASER EQUATIONS 

In the previous section we have derived the laser equa- 
tions in the uniform field limit for a resonant laser. For 
arbitrary detuning it can be demonstrated that the laser 
equations in the uniform-field limit read |l6j | 



dF 
~dt 



dF 

oz 



k(P-F) 



dP 
~dt 



: 7X [FD - (1 + iA c ) P] 



dD 
~dt 



711 [r - D - \ (FP* + F*P)] 



(157) 
(158) 
(159) 



where A c = (uj c — LU21) /7_l is the atom-cavity detuning 
parameter. These equations are complemented with the 
periodic boundary condition 



F (0,t) = F (L m ,t) 



(160) 



The fact that the boundary condition is periodic means 
that the field can be written in the form 



F(z,t) 



+ 00 

E 

m— — oo 



F m (t) , 



(161) 



where q m — m27r/L m with m an integer. This means 
that the intracavity field is, in general, a superposition of 
longitudinal modes of the empty (i.e., without amplifying 
medium) cavity. In order to see this clearly, let us solve 
Eq. (|157H for the empty cavity and ignoring cavity losses, 
i.e., Eq. (|157fl with its right-hand side equal to zero. Its 
solutions have the form 



F(z,t) 



+ 00 

E 

m— — oo 



(162) 



with tu m = vq m . Now we must notice that the actual 
field is not F (z, t) but F (z,t — t) with r = zAt/L m and 
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At = [L c — L m ) /c, as we introduced new fields in Eqs. 
I(126[l and 1127(1 (remind that the field F is proportional 
to the field a' which is different from the actual field 
a). Then, the actual field (we will not introduce a new 
symbol for it) is 



read 



+ oo 



F(z,t)= F ™ ex P 



L c — L, 



i q m z - ijj m z- 



+00 



= ^ F m exp [i (k m z - uj m t)} , 



ui m t 

(163) 

(164) 



with 



w m L c — L m 2ttc ,.„., 
k m = q m p = m— , (165) 

C L m Lie 



where ui m = vq m has been used. Eq. I|165|) shows clearly 
that the actual field appears decomposed into empty- 
cavity modes. 

Thus Eqs. l|157H159|l can describe multilongitudi- 
nal mode emission when dF/dz is non null and, when 
dF/dz = 0, this model can describe only singlemode 
emission. The question now is: Should we keep the spa- 
tial derivative always? Or, in other words, when will the 
laser emit in a single mode and when in several longitudi- 
nal modes? In 1968 Risken and Nummedal and. in- 
dependently, Graham and Haken [181 ] , demonstrated that 
Eqs. l|157H159]l predict the existence of multilongitudi- 
nal mode emission if certain conditions are verified. We 
are not going to treat the Risken-Nummedal-Graham- 
Haken instability here (see, e.g., 0, H E3 13 or [3 
for a recent review), it will suffice to say that for mul- 
tilongitudinal mode emission to occur the two necessary 
conditions are: (i), a large enough pump value (in reso- 
nance, A = 0, r must be larger than nine and remember 
that the laser threshold in these conditions, given Eq. 
(|105ll . equals unity; out of resonance even more pump is 
required), and most importantly; (ii), the cavity length 
must be large, unrealistically large for common lasers. 
Then for short enough cavities (and this is not a re- 
strictive condition at all for most lasers) the laser will 
emit in a single longitudinal mode. In this case, the spa- 
tial derivative in Eq. I|157f) can be removed and we are 
left with the Maxwell Bloch equations for a singlemode 
laser. We must insist that all this is true for homoge- 
neously broadened lasers and cannot be applied to inho- 
mogeneously broadened ones, see [l6| . 

Then for singlemode lasers we can take dF/dz = 0. It 
is particularly interesting to write down the singlemode 
laser equations in resonance (A c = 0). Let us write the 
field and atomic polarization in the following way 

F = Ee**, P = {Pre + iPim) e** , (166) 
with E a real quantity. Now Eqs. I|157I158|) . with A c = 0, 



E = k {P re - E) , 
p. 

, * irn 

4> = k— , 

Pre = 7± (ED -Pre)+ 
P%m — Ti-^m 



(167) 

(168) 

(169) 
(170) 



where the dot means total derivative with respect to time. 
By suitably combining the first, the second and last equa- 
tions one obtains 



P E 

p- + p = -(K + 7J.) 



from which 

Pirn \p) 



Pirn (0) 
25(0) 



E{t)e 



-(K+7J.)t 



0. 



(171) 



(172) 



and then (j> — * also. Thus, in resonance, the singlemode 
laser equations reduce to only three real equations 



E = k(P-E), 
P = 7 _l (ED - P) , 
D = 7|| (r- D - EP) 



(173) 
(174) 
(175) 



where P = P re - 

The above set of equations is usually known as Haken- 
Lorenz equations. The reason for this name is the follow- 
ing: Let us define the adimensional time r = 7_i_i, and 
the new variables and normalized relaxation rates 



X = E, Y = P, Z = r-D, 

k 7|| 
a = — , = — . 

7-L 7-L 



These new variables verify 



-X = a(Y-X), 

—Y = rX -Y + XZ, 
dr 

— Z = b(Z - XY) . 
dr 



(176) 
(177) 



(178) 
(179) 
(180) 



These are the Lorenz equations |l9j , which are a very sim- 
plified model proposed by Edward N. Lorenz in 1961 for 
the baroclinic instability, a very schematic model for the 
atmosphere. They constitute a paradigm for the study 
of deterministic chaos as they constitute the first model 
that was found, by Lorenz himself, to exhibit determin- 
istic chaos. It was Herman Haken who, in 1975 
demonstrated the astonishing isomorphism existing be- 
tween the Lorenz model and the resonant laser model 
that we have just demonstrated. After this recognition, 
the study of deterministic chaos in lasers became a very 
active area of research (see, e.g., 0, 0,11 El 0)- 
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The Haken-Lorenz model exhibits periodic and chaotic 
solutions, and several routes to chaos can be found in its 
dynamics. The equations can be numerically integrated 
easily, e.g. with Mathematica, and we refer the inter- 
est reader to |2j, LD, |9j, |l£J, [l2j for suitable introductions 
into these fascinating subjects. Here we shall only briefly 
comment on a particular point. 

Eqs. (|173I175[) have two sets of stationary solutions: 
The laser off solution (E = P = and D = r), and the 
lasing solution (E = P = ±\Jr — 1 and D = 1) that 
exists for r > 1. A linear stability analysis of this last 
solution shows that it becomes unstable when n > "f± +71 
(this condition is know as "bad cavity" condition) and 
r > thb with 



TUB 



k(k + 37j_+7||) 

K - (7_L +7||) 



(181) 



whose minimum value is 9 for 7|| = and k = 
Notice that as we are considering the resonant case for 
which the lasing threshold is r on = 1, Tub — 9 means 
that the adimensional effective pump r must be, at least, 
nine time above the instability threshold. 

From Eq. (|181fl we can say that the singlcmode so- 
lution is always stable for good cavities (k < 7x + 7||) 
and also for bad cavities if the pump is small (r < Thb), 
but for bad cavities and large pump, the stationary solu- 
tion becomes unstable (through a Hopf bifurcation) and a 
self-pulsing occurs (e.g., chaotic oscillations). The con- 
dition r > tub (remember that Thb > 9) is usually 
considered as a very restrictive condition (we insist, the 
laser should be pumped nine times above threshold, and 
this is quite a large pump value!) but we have seen that 
we must be careful when interpreting the pump param- 
eter r. In fact, if one considers a three-level laser and 
uses Eqs. <|103|) . i|H2[) . (G>SJ) an d (EH, one can write Eq. 
(|181fl in terms of the actual pump strength and decay rate 
and gets that the instability threshold to lasing threshold 
pumps ratio can be very close to unity. 

We can show this easily. Let us recall equation (|l(J3f) . 
that relates the adimensional effective pump parameter 
r with the inversion in the absence of fields do 



r = Gd 

2gL n 



G = 



7_Lc|ln^T 



(182) 
(183) 



as well as Eq. I|54[) that relates do to the actual physical 
pump parameter R in three-level lasers. By taking into 
account that r on = 1 and tub = 9 it is easy to see that 



R 3 h L b _ (G + 9)(G-1) 
Rlk (G + l)(G-9)' 

which for large G simplifies to 



(184) 



(185) 



i.e., for three-level lasers the "very restrictive" condition 
Tub = 9 turns out to be an easy condition in terms of 
pumping {Rife/Ron = 1 + e) when the gain parameter 
G, Eq. (|183fl . is large enough. 



IX. CONCLUSION 

In this article we have presented a self-contained 
derivation of the semiclassical laser equations. We have 
paid particular attention to: (i) the adequacy of the stan- 
dard two-level model to more realistic three- and four- 
level systems; and (ii), the derivation of the laser equa- 
tions in the uniform field limit. We think that our pre- 
sentation could be useful for a relatively rapid, as well as 
reasonably rigorous, introduction of the standard laser 
theory. This should be complemented with a detailed 
analysis of the stability of the stationary laser solution 
(see, e.g., UGHUElEl]) which we do not treat here for 
the sake of brevity. 
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X. APPENDIX 

In this Appendix we demonstrate that the adiabatic 
elimination of the atomic coherence in Eqs. I|41I42() con- 
sists in making dt<J\2 = 0. 

Consider the evolution equation 



^/(*) = -7/ (*) + <?(*), 



(186) 



which must be complemented with the equation of evo- 
lution of g (t) . Notice that Eq. (|186[) coincides with Eq. 
(021 f° r f (t) = c\i exp {iSt), g{t) = ida* exp (i5t) , and 
7 = 7^. Now we define the new variable f — f exp (—"ft) 
that verifies 



dt 



f 



-ft 



9®, 



from which 

/(*) = /(0)+ / ,//', -'>/(/'>. 

and then 

/(«)=/ (0) e-* - 



dt'e-^-^git') 



(187) 



(188) 



(189) 



Integrating by parts and ignoring the first (decaying) 
term 



/(*) = - 

7 



git) 



dt'e-<^')±g(t') 



(190) 
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Then, after repeatedly integrating by parts, one finally 
obtains 



7 1 g(i), which is the result one obtains by making 
4zf (t) = in Eq. (|186[) . as we wanted to demonstrate. 



Id Id 2 

1 ~ 7d7 + y di 2 " 



g(t), (isi) 

and thus for large enough 7 one can approximate / (t) « 
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